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Abstract 

We study the dynamics of electroweak symmetry-breaking in an extension of the Standard Model where the Higgs sector is augmented by 
the addition of a real (Y = 0) isospin triplet. We show that this scenario exhibits a novel, two-step electroweak phase transition, wherein the 
first step provides the strongly first order transition as required for electroweak baryogenesis followed by a second step to the Standard Model 
Higgs phase that also admits a cold dark matter candidate. We analyze the constraints on this scenario from recent results obtained at the Large 
Hadron Collider for the Higgs diphoton decay channel. We argue that this two-step scenario can be generalized to extensions of the Standard 
Model with additional higher-dimensional scalar multiplets that may yield realistic baryogenesis dynamics. 



I. INTRODUCTION 

An outstanding problem at the interface of cosmology with 
high-energy and nuclear physics is to explain the origin of 
baryon asymmetry of the universe: 

Fb = — = (8.82 ± 0.23) X 10"" WMApH (1) 
s 

where ub and s are the baryon and entropy densities, respec- 
tively, and where the value has been determined from studies 
of the cosmic microwave background. Assuming a matter- 
antimatter symmetric initial conditions, three ingredients must 
have been present in the microphysics of the early universe to 
generate a non-vanishing Yb 121: (1) violation of baryon num- 
ber; (2) violation of both C- and CP-symmetry; and (3) depar- 
ture from equilibrium dynamics or CPT violation. An attrac- 
tive mechanism for solving this problem is electroweak baryo- 
genesis (EWB) (for a recent review and extensive references, 
see Ref. lUtl), wherein baryon number generation is driven by 
the generation of CP asymmetry at the time of the electroweak 
phase transition (EWPT). An important ingredient for its suc- 
cess is the existence of a strong first order electroweak phase 
transition in which electroweak symmetry-breaking (EWSB) 
proceeds via bubble nucleation. 

It is well known that the scalar sector of the minimal Stan- 
dard Model (SM) is unable to have generated a strong first 
order phase transition because the Higgs boson is too heavy'. 
Extensions of the scalar sector of the Standard Model have 
been motivated, in part, by alleviating this difficulty. Ap- 
proaches include (a) the introduction of new scalars that in- 
crease the barrier between the broken and unbroken phases 
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associated with the non-analytic 0{T) in the finite tempera- 
ture effective potential; (b) the introduction of tree-level cubic 
terms in the T = potential that yield a barrier between the 
two phases; and (c) reliance on logarithmic corrections to the 
potential that can drive the transition (for a recent discussion, 
see Ref. |4J). 

Most analyses of the EWPT have thus far relied on only 
one of these approaches while making the reasonable - but 
not necessary - assumption that EWSB proceeds in a single 
step. It is possible, however, that the dynamics of EWSB are 
more complicated, involving more than one of the aforemen- 
tioned mechanisms and proceeding through multiple steps, 
passing through intermediary phases before reaching the EW 
phase. In this work, we study a minimal scalar sector ex- 
tension that gives rise to a two-step transition and that re- 
lies on two of the general mechanisms listed above. Specif- 
ically, we add a real triplet E that transforms as (1, 3, 0) un- 
der SU(3)c X SU(2)l X U(I)y, corresponding to the small- 
est dimension scalar multiplet that carries non-trivial SU(2)l 
charge. Using this minimal scenario, we show how a multi- 
ple step transition may facilitate the first order EWPT needed 
for EWB, yield a dark matter particle, and contain dynam- 
ics testable through Higgs boson decays into diphotons at the 
Large Hadron Collider 

In brief, at temperature T ~ the neutral component T,^ 
can provide a viable cold dark matter candidate when the 
scalar potential admits a Z2 symmetry (for an analysis of the 
zero-temperature properties and collider phenomenology, see 
Ref. iH). At r 100 GeV, the interplay of the scalar triplet 
and Higgs doublet fields (in the minimum of the free energy) 
may give rise to a two-step transition, schematically illustrated 
in Fig. [T] During the first step, at a temperature just below 
the first critical temperature T^, the system makes a transition 
from the symmetric phase at point O to the isospin breaking 
phase at point S, where the neutral component of the triplet 
field obtains a non-vanishing vacuum expectation value (vev) 
while the neutral Higgs vev remains zero. Then, in the 
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second step at a temperature below the second critical temper- 
ature Th, the universe makes a transition from isospin break- 
ing phase at point S to the electroweak symmetry breaking 
phase at point H, where the Higgs vev {H'^) becomes non- 
vanishing but (S") relaxes to zero, ultimately leading to the 
T = Higgs phase with as the dark matter particle. 

While the remainder of the paper addresses the dynamics 
for this scenario in detail, we comment on several salient fea- 
tures here. 

(1) Although the zero temperature I]° vev in the Higgs 
phase need not vanish, the constraint from the elec- 
troweak p parameter requires it to be small < 
4 GeV. Allowing it to be tiny but non-vanishing does 
not substantially alter the EWPT dynamics but does 
preclude as a viable dark matter candidate. Con- 
sequently, we take (S") = at the conclusion of the 
second step in order to yield a dark matter candidate. 
To that end, we impose a dark matter-preserving (^2)2 
symmetry on the potential and refer to the model as the 
"Z2SSM." 

(2) The first step of the phase transition (O — > E) can be 
strongly first order, driven entirely by the finite- T dy- 
namics of the effective potential along the neutral triplet 
scalar direction. This transition is analogous to the one 
that might occur along the O ^ H direction in the 
Standard Model but is excluded from being first order 
due to non-existence of a sufficiently light Higgs boson. 
In contrast, the parameters in the S sector of the theory 
are sufficiently unconstrained by current phenomenol- 
ogy to allow for a strong first order EWPT along the 
O ^ E direction. 

(3) Bubble nucleation during the first step (O — > S) creates 
the necessary environment for baryon number genera- 
tion, assuming additional sources of CP-violation be- 
yond those provided by the Standard Model. Since 
the E-field carries non-trivial SU(2) l charge, the Bh-L- 
violating sphaleron transitions that create and destroy 
baryon number are quenched inside the S-phase bub- 
bles, capturing any net baryon number density produced 
ahead of the advancing bubble walls. In the present 
work, we concentrate on the phase transition dynam- 
ics for this step, leaving an analysis of possible sources 
of CP-violation to a future study. 

(4) The dynamics of the second step (S H) are gov- 
erned by the tree-level interaction between the S and 
H fields. To ensure that any baryon number produced 
during the first step is not washed out by reactivation 
of the electroweak sphalerons, this transition must also 
be both first order and sufficiently strong. Moreover, 
entropy production cannot be too copious without di- 
luting the initial baryon asymmetry^. As we discuss be- 




FIG. 1 : Field phase space indicating critical (extremal) points in the 
tree-level potential at zero temperature, and the expected two-step 
pattern of symmetry breaking at finite temperature. Red points are 
related to black points via {Z2)h and {Z2)t. symmetries. 



low, the degree to which these requirements are met can 
be constrained by measurements of the 77 rate. 

(5) While the introduction of the E field constitutes the 
minimal extension of the Standard Model scalar sec- 
tor leading to this scenario, it is possible that higher- 
dimensional scalar multiplets will yield similar dynam- 
ics, though not necessarily with a cold dark matter can- 
didate as well. As we argue below, the main features 
are otherwise generic and are likely to persist in other 
models that also include additional degrees of freedom 
as needed for an appropriate ultraviolet completion. In 
this regard, a similar two-step scenario was considered^ 
in |j6t| involving a second Higgs doublet rather than a 
real triplet. The authors envisioned baryon-asymmetry 
generation to occur during the much stronger second 
step of the phase transition. However, as pointed out in 
101, since weak isospin symmetry is already broken in 
the first step, Bh-L violating processes are already too 
suppressed to convert existing CP asymmetry to baryon 
asymmetry. Hence, we concentrate on the case where 
the first step is strongly first order 

Our analysis of the foregoing scenario is organized as 
follows: in section [III we formulate the model and subse- 
quently discuss the zero temperature vacuum structure and 
tree-level vacuum stability constraints in section |III] We then 
turn our attention to the finite-T dynamics, focusing first on 
the sphaleron transitions in section |IV] then following their 
behavior during the two-step EWPT in section |V] Finally, we 
draw connections to collider phenomenology and discuss im- 
plications coming from recent LHC results in section |VT] 
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II. MODEL AND FORMULATION 



III. ZERO TEMPERATURE VACUUM STRUCTURE 



In the real triplet extension of the Standard Model, the 
scalar sector is composed of the standard Higgs doublet H ^ 
(1,2, +1/2), and a real triplet^ E - (1,3,0), with Q = 
+ Y . Our convention for component fields are 



H 




where we fix the tree-level Higgs vacuum expectation value 
(VEV) Wo = 246 GeV and where xq = (S") is the triplet vev. 
We denote the physical quanta of charged and neutral fields 



= {a I =F i(T2)/v2 and — a^, respectively. 
Following the notation of our earlier workljH we write the 
Z2-symmetric scalar potential as 



V{H, S) = -fi'^H^H + \{H^Hf - ^(S • E) 



64 



4 ' ' 2 



(S-S)^ 



WH{Y. ■ E) 



(3) 



This potential exhibits two independent global symmetries 
that are also accompanied by reflection symmetries^: 



SO(4)h : H ^ Oi^^H 
SO(3)s : S ^ O3X3S 



{Z2)h : H^-H 
(Z2)s : -E. 



We use these symmetries to restrict the discussion of the elec- 
troweak phase transition to electrically neutral components of 
the Higgs doublet and isospin triplet. We should mention, 
however, that as detailed in f5J, the addition of a possible 
ai7J^(E • T)H breaks the SO(3)s symmetry and would trig- 
ger the acquisition of a vacuum expectation value for the neu- 
tral component of the isospin triplet, restricted to values of 
the order (E") < 4 GeV < (H) = 246 GeV by the tightly 
constrained electroweak p-parameter at tree-level. Although 
its presence may substantially affect collider phenomenology 
even when within bounds of experimental uncertainty for pEW 
through its mixing with the doublet, the impact on the elec- 
troweak phase transition is negligible. Therefore, for clarity, 
we restrict our discussion of the phase transition to the region 
of the parameter space where the (^2)1; symmetry is realized, 
thereby also allowing the neutral component to be a dark mat- 
ter candidate. 

This model contains only three real parameters in addition 
to those of the Standard Model: the triplet (negative) mass- 
parameter triplet self-coupling 64, and the doublet-triplet 
"Higgs-portal" coupling 02. The analysis is greatly simpli- 
fied due to the addition of small number of undetermined pa- 
rameters and the limited number of field degrees of freedom 
participating in the electroweak phase transition. 



We follow the convention that the arrow above a letter denotes a vector in 
isospin space. 
' The Yukawa sector of the Standai 
symmetry for the Higgs doublet. 



isospin space. 

^ The Yukawa sector of the Standard Model Lagrangian breaks the reflection 



The qualitative behavior of the electroweak phase transition 
is largely influenced by the zero temperature vacuum structure 
of the Higgs sector In this section we perform a vacuum sta- 
bility analysis at tree level. Following our conventions in (|2]i, 
we denote the neutral components of the isospin doublet and 
triplet as h and a respectively. Upon setting the remainder of 
the components to zero, the potential reads 



1 



^ 2 2 , 1 \ 1,4 



(4) 



The critical points are found by solving the minimization con- 
ditions 



dVr, 



dh 



9K< 



da 



0. 



(5) 



The potential, being a fourth-order polynomial in two vari- 
ables, nominally has a total of nine extremal points. Apart 
from the origin (/i, a) = (0, 0), the extrema come in posi- 
tive/negative pairs related by combinations of the {Z2)h and 
(■^2)2 reflection symmetries. We use these symmetries to 
eliminate the redundant negative partners down to four dis- 
tinct critical points as indicated by black points in Fig. [T] la- 
beled by H, O, E, and T. We note that not all points may 
be realized as critical points of the potential as certain choices 
of parameters may yield complex-valued solutions (see poten- 
tials in Fig. O. 

We label point H as the Higgs phase point that corresponds 
to the physical electroweak vacuum, where the tree-level mass 
relations are determined to be 



to|± 



(125 GtW f (LHC) 
1 9 



(6) 



More generally, the requirements of vacuum stability may ef- 
fectively be summarized by the condition that point H, the 
electroweak vacuum located at {h, a) = (246 GeV, 0), is the 
global minimum of the potential, with vri^ > and ni^^ > 
masses. This requirement is expressed by the inequality 



1 /'I 2 2\ 1 



(7) 



Further, in order to facilitate the discussion of two-step 
phase transitions, it will be useful to identify regions of param- 
eter space where the potential exhibits a secondary local min- 
imum at point E with positive masses. A straightforward cal- 
culation yields the condition for the existence of a secondary 
minimum. 



1 2 1 a2 / 1 2 



(8) 



which requires /i| > in (|6]l. 

In Fig. |2]we display the region (shaded yellow and blue) 
in the 02-^4 plane for which the vacuum stability condition in 
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(|7]l is satisfied, with the masses tos = 150 GeV, and niH = 
125 GeV held fixed. The blue shaded region indicates points 
where the requirement (|8} is also satisfied and the potential 
has a secondary local minimum at point S. To assist the reader 
in visualizing the potential for various regions of parameter 
space, we provide illustrative plots in Fig. |3]of the potential 
for two cases: (a) Eq. (Q only being satisfied corresponding 
to a representative point in the yellow region in Fig. |2]and 
(b) both Eqs. (|7]l and ^ holding, corresponding to the blue 
region in Fig. |2l 



A. Quantum corrections 

Quantum corrections to the tree-level potential can, in gen- 
eral, have a significant impact on the vacuum structure of the 
theory. To obtain a reasonably reliable analysis of loop ef- 
fects, it is useful to employ a renormalization group (RG) im- 
proved potential, wherein one replaces the parameters by the 
corresponding running parameters that depend on the renor- 
malization scale /i. Choosing /i commensurate with the value 
of the field effectively re-sums large logarithms that would 
otherwise impair the perturbative expansion, moving the log- 
arithms into the running parameters and yielding explicit cor- 
rections that are manifestly perturbative. In general, the use 
of running parameters can lead to either instabilities or non- 
perturbative couplings at high scales, depending on the con- 
tributions to the /3-functions. In the SM, for example, the top 
quark loop contributions may cause Higgs quartic self cou- 
pling to run negative below the Planck scale if the input value 
at the electroweak scale - set by the value of the Higgs mass 
- is too small (for recent work on the implications for vacuum 
stability and perturb ativity in scenarios for physics beyond the 
SM, see e.g. Refs. lHH)- 

In the present instance, the potential remains stable along 
the cr-direction up to the Planck scale since the vanishing 
hypercharge of the S precludes any couplings to fermions 
and, thus, any substantial negative contributions to the /3 
function for the E quartic self-coupling, 64. Moreover, the 
WHF interaction in ^ generate a positive contribution to 
px proportional to a^, partially compensating for the negative 
top quark contribution and improving stability of the Higgs 
and S phases. On the other hand, A(/i) may become non- 
perturbative below the Planck scale if 02 is too large. Since 
our interest here focuses on the novel phase transition dynam- 
ics associated with this scalar sector extension, we defer to 
future work an analysis of the scale at which the Z2SSM be- 
comes non-perturbative. 



generating step. Within the symmetry-broken bubble, a vari- 
ant of the Standard Model sphaleron, which we term the E- 
Sphaleron, is active and threatens to erase any baryon asym- 
metry. In this section, we analyze the the S-sphaleron con- 
figuration that is specific to the model, compute its rate in a 
gauge-independent manner at finite temperature to estimate 
the condition on the strength of the transition needed to avoid 
baryon number washout. 

The degree to which an initial baryon asymmetry is erased 
depends on the sphaleron transition rate inside the SU(2)i- 
broken phase: 



r,ph = A{T) cxp i-Esph/T) 



(9) 



where E'sph is the energy of the sphaleron configuration 
and A{T) depends on the frequency associated with the 
sphaleron unstable direction, zero modes associated with ro- 
tational and translational symmetries, and fluctuations about 
the sphaleron^. Roughly speaking, the fraction of the initial 
baryon asymmetry that remains after the EWPT completes 
goes as exp(— Fsph). Thus, the larger the sphaleron energy, 
the slower the washout rate and the more effectively baryon 
number is preserved. A detailed analysis of baryon num- 
ber washout then leads to a "baryon number preservation cri- 
terion" (BNPC) that, in turn, depends on a computation of 
the s pha leron energy (for a recent, detailed discussion, see 
Ref. iJitl). In the SM, the BNPC can be expressed as a re- 
quirement on the ratio of a finite temperature scale v{T) (de- 
fined below) to T at the critical temperature for the transition. 
To determine the corresponding BNPC for the Z2ESM, we 
require a computation of the sphaleron energies for the two 
steps of the transition. 

To that end, we start with the zero-temperature energy func- 
tional for the EW sphaleron, given by 



E = 



where 



-W,j ■ W, 



^3 



—BiiB 



{d,h)Hd,h) 



-Tr[(AS)'] + y(i?,S) , 



(10) 



w, 



d^Wj - djWt - gW^ X Wj 
Bij = diBj — djBi 



DiH = {d, + igf ■ + ig'YB,)H 
D,t = d^t + .gS X Wi 

and V{H, S) is the tree-level potential of (|3]l. The generaliza- 
tion to finite temperature is treated in later in this section. 



IV. E-SPHALERON AND BNPC 



A. Ansatz 



Crucial to baryogenesis in this model is the two step pro- 
cess by which the phase transition proceeds. In the first step 
the system makes a transition from the symmetric phase to 
the weak-isospin broken S phase at a critical temperature 
we define as T^,. This step is to proceed via bubble nucle- 
ation since we expect this step to be the baryon-asymmetry 



Our ansatz at spatial infinity for the isospin triplet meson 
is motivated by the Klinkhamer-Manton lflTll ansatz applied to 



' Assuming -Egph is computed in the saddle point approximation. 
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m«=125 GeV, ^2=1.07 m„=l25 GeV, a2=\.01 




100 120 140 160 180 200 -2 -1 1 2 

wis 

FIG. 2: Regions A (yellow striped) plus B (solid blue) indicate where tree-level electroweak vacuum stability condition of l|7j is satisfied. Left 
panel: the ms-bi plane for fixed (mn — 125 GeV, a2 = 1.07); right panel: the 02-64 plane for fixed (m^r = 150 GeV, jns = 150 GeV). 
Regions B indicate where l[8j is also satisfied and the tree-level potential exhibits a metastable minimum along neutral E direction. Illustrative 
representations of the scalar potential for regions A and B are indicated in the left and right panels of Fig. |3] respectively. 




FIG. 3: Qualitative picture of the potential V{h, a) of ^ in the two different regions of parameter space as indicated in Fig. (2] Potential 
A (corresponding to regions A of Fig. |2) displays no critical point along the a direction, whereas Potential B (corresponding to regions B of 
Fig.[2l( exhibits a metastable minimum along the cr-direction. 



the isospin doublet, in which the field configuration at spa- 
tial infinity may be viewed as a pure gauge transform of the 
vacuum of the form 

U{x) = Xia.i, (11) 

so that at spatial infinity the isospin triplet scalar S(x) = 
S(x) • T and gauge fields W(x) = W(x) • T approach 

I](|x| ^ CX)) C/(x)SvacC/^(i) 

= (X • S]vac)(i ■ <?) - ^Svac ■ (? (12) 



W(|x| ^oo) = =lU{x)\JU'^{x) 



9|x| 



X X a . 



(13) 



where Eyac — (0, 0, x-q) is taken to be the homogenous vac- 
uum configuration minimizing the potential V{H, E) at the 
metastable point S in Fig. [T]along the triplet direction. Thus, 
our ansatz for the isospin-triplet S(x) and gauge field W(x) 
at all spatial points is given by 



E(x) = s(x) [{x ■ Svac)(^ • ct) - -Evac ' o"] 

W(x) = X a 

9m 



(14) 



where s(x) and /(x) stand for the sphaleron radial profile 
functions to be determined by minimizing (fTOl i subject to the 
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boundary condition 



s(0) = 0. 
/(0)-0, 



s{oo) 

/M- 



(15) 



A substitution of ( fT4b into the energy functional ( fTOl i yields 
the functional dependence on radial profile functions /(f) and 
s(f): 



sph 



47rxo 



92 Jo 



dr 



2 



1 6, 



3 



Attxq 

92 



BT,SM{b4/92 



(16) 
(17) 



where we have defined f = (722;o|x| the dimensionless ra- 
dial coordinate measured from the center of the sphaleron. We 
note the functional dependence on the radial functions is iden- 
tical to that of the Klinkhamer-Manton sphaleron, apart from 
the numerical coefficient of the coupling term s^(l — /)^ in- 
volving the gauge and scalar field profile functions being 8/3 
instead of 1. 



B. Numerical minimization 

In this section, we numerically minimize the energy func- 
tional ( fT6] l crudely following the variational principle. We 
make use of the following ansatz for the radial profile func- 
tions: 



fir) = r^— 2 
s(f) 



(18) 



where the characteristic sizes pj and stand in for the vari- 
ational parameters. 

In Table m we summarize the results of the numerical min- 
imization of the energy function Bj^sm for select values of 
the coupling ratio 64/.g|. For comparison, also displayed are 
the results for the Standard Model sphaleron energy function 
BsM for corresponding coupling ratio A/gl using an identical 
ansatz to that of (fTSl l . The difference in energies of the B 
function of the two sphalerons is entirely represented by the 
differing coefficients of the coupling term .s^(l — /)^ in ( fT6] l. 



C. Baryon number preservation bound 

So far, the discussion of the sphaleron energy has been re- 
stricted to zero temperature. To obtain the sphaleron rate at 
finite temperature, the tree-level potential V{H, S) appearing 
in the energy functional (fTOl i is replaced with the temperature 
dependent effective potential V{H, E; T) in the high-T ap- 
proximation to ©(T^), shown in dZST l below. The result is that 



Sigma-SM sphaleron 


Standard Model sphaleron 




Pf 


Ps 


5esm 




Pf 


Ph 


BsM 


0.001 


2.274 


5.273 


2.015 


0.001 


2.862 


4.120 


1.619 


0.01 


1.999 


3.714 


2.156 


0.01 


2.609 


3.058 


1.693 


0.1 


1.654 


2.101 


2.485 


0.1 


2.250 


1.752 


1.883 


1. 


1.373 


1.057 


2.969 


1. 


1.962 


0.856 


2.155 



TABLE I: Results of the numerical minimization of l il6b using the 
ansatz J18t for specific values of ratios of coupling constants for the 
left: E-sphaleron and right: Standard Model sphaleron. For a given 
ratio of the coupling constants, the E-sphaleron represents a greater 
barrier than that of the Standard Model (see discussion in text), lead- 
ing to a relative suppression of the sphaleron rate. 



the sphaleron energy at finite temperature is scaled Hi 811 rela- 
tive to the zero temperature energy; 



i?sph(r) - ^i?sph(o) 



x{T) ( AttB 



SSM 



2^0 V 92 



(19) 



where x{T) is the gauge-independent finite tempera- 
ture expectation value (S") derived from minimizing the 
V{H, S;r), and -Esph(O) is the zero-temperature sphaleron 
energy numerically computed above. A similar expression ap- 
plies to the energy of the SM sphaleron. 

The actual strength of the phase transition necessary to pre- 
serve the baryon asymmetry requires a calculation for the gen- 
eration of baryon asymmetry is to some extent model lll6ll . 
Both an actual calculation of the magnitude of the generated 
baryon asymmetry and an estimate of the semi-classical pref- 
actor of the sphaleron rate A{T) are required for placing a 
requirement for the phase transition strength. Since neither 
of these calculations are presently available, we make a rough 
estimate of the bound by comparing the numerical results of 
corresponding B functions that enter into the sphaleron en- 
ergy. For 0(1) values of the coupling ratio h^/ g^, the ratio of 
the B functions (from Table [III is 



gsSM(l.O) 

Ssm(I.O) 



4 
3 



(20) 



which we expect to be fairly insensitive to the ansatz used in 
(HI. 

From (fT9] l we then expect that the BNPC bound on the 
"order parameter" of the first step of the phase transition 
(O — > S) is relaxed approximately 25% relative to that of 
the Standard Model bound: 



first step BNPC: 



> 0.75, 



(21) 



A similar bound on v{T)/T should be placed on the second 
step of the phase transition (S H), where v{T) is the finite 
temperature Higgs vev. We expect the bound to be similar that 
in the Standard Model 



second step BNPC: 



T^nuc 



> 1.0, 



(22) 



but is subject to be refined by calculations of the baryon asym- 
metry and sphaleron rate prefactor For the analysis of the 
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phase transition in the next section, we shall loosely take the 
BNPCs (I2TI) and (|22]) to guide us in parameter space to where 
a successful electroweak baryogenesis scenario could be real- 
ized. 



V. EW PHASE TRANSITION 

We now follow the dynamics of the EWPT in the Z2i;SM. 
Before doing so, we comment on the issue of gauge invari- 
ance that has been at topic of recent interest in this context. 
One of the unresolved theoretical issues plaguing perturbative 
analyses of the electroweak phase transition is the problem of 
extracting gauge-invariant quantities relevant to baryogenesis. 
The root of the problem lies in the lack of a gauge-invariant 
definition of the free energy that is compatible with pertur- 
bation theory. Although the problem has been known since 
the early days of the development of thermal gauge theories 
lfl9il20ll . it is in our view that the problem not yet been solved 
satisfactorily. The choice of gauge most commonly employed 
in the literature is the 't Hooft background Landau ^ = 
gauge. 

In an earlier paper by us III6I1 . the issue was tackled, and 
a theoretically sound method to carry out an analysis that 
is fully gauge-independent and compatible with perturbation 
theory was developed. While the method applied to the SM 
as well as some features of the minimal supersymmetric SM 
gives the correct qualitative dependence on model parame- 
ters when compared against lattice results, maintaining gauge- 
independence at the perturbative level appears to come at the 
expense of numerical accuracy. 

With this caveat in mind, the attitude we adopt here is the 
moral one: to compute physical observables in a way that suf- 
fer from no gauge ambiguities, and so we will follow the pro- 
cedure detailed in our earlier paper with the awareness of the 
method giving correct parametric dependence but not numer- 
ical accuracy. We now proceed using the gauge-independent 
procedure as outlined in [fl6ll for finite temperature analysis of 
the free energy of the system. 

The construction of the finite temperature one-loop effec- 
tive potential follows the standard background-field method 
ll2lll . We write the resulting effective potential schematically 
as: 

V{h,a,T) = Vn-,,{h,a) 

+ h{V^=\h,cj) + V^*\h,a,T)) , (23) 

where Vaee{h,a) is the tree-level potential in (|4]l, and 
V^^'^{h,a) and Vi^^{h,(7,T) are the zero-temperature 
Coleman- Weinberg and temperature-dependent potential. Ex- 
plicit expressions in the 't Hooft background-field gauge 
is given in appendixlAl 

We follow the evolution of the minima of the effective po- 
tential as a function of temperature by inserting the critical 
values of the tree-level potential defined by (|5]l into the one- 
loop temperature-dependent potential (|23] l. at which point the 
gauge dependence cancels. As already discussed above, there 
are generally four such critical values to follow, corresponding 



to phases at H, E, O, and T in Figure[T| We determine critical 
temperatures by requiring that the degeneracy condition 

= (24) 

be met. 

In Fig. m we provide two representative examples for the 
evolution of the free energy at its critical points as a function 
of temperature. In Fig. |4^, we choose a large triplet scalar 
mass ms = 170 GeV, corresponding to a point inside the 
yellow region of Fig. |2^, where the zero-temperature potential 
takes the form as in Fig. [3^ with no critical point along the a 
direction. The evolution of the phase structure with tempera- 
ture corresponds to a trajectory that follows the curve of low- 
est energy. Thus, at high-T the universe proceeds along the 
(blue) line labeled O, corresponding to the symmetric phase. 
For T <Th ^ 108.4 GeV, the (red) cui-ve labeled H, giving 
the value of the potential at the electroweak minimum, has the 
lowest energy, so the trajectory then switches to the (red) H 
curve. A first order transition to the H phase will then occur 
at T just below if the nucleation probability is sufficiently 
large. In the SM, this transition occurs via cross over rather 
than bubble nucleation. 

On the other hand in Fig. [3]5, for a lighter triplet mass 
jTis = 130 GeV, corresponding to a point in the blue region 
of Fig. [2^, the zero-temperature potential takes the form as 
in Fig. [3}5, where a metastable minimum along the a direc- 
tion exists. Consequently, the system exhibits a richer two- 
step phase transition, favorable for baryogenesis, with critical 
temperatures at T„ = 122.7 GeV and T„ = 86.7 GeV. The 
quantity relevant to the preservation of the baryon asymmetry 
at the time of the electroweak phase transition is the order pa- 
rameter at the nucleation temperature: 2;(r™'^)/T'""'^ for the 
first step, and for the second step, as defined in 

(ISTT i. We use the more conservative order parameter, evalu- 
ated at the critical temperature, lying somewhat higher than 
the nucleation temperature. 

In Fig. |5] we display the results of a scan for fixed 02 = 
1.07 and mn ~ 125 GeV, which is superimposed over the 
tree-level vacuum stability graph of Fig. |2] The contours are 
shown only where a two-step phase transition occurs. The 
red contours correspond to lines of constant order parameter 
x{Tcr)/Tcr, relevant for the first step, and the blue contours are 
lines of constant v{Th)/Th for the second step. We observe 
that the strength of the transition in the first step, character- 
ized by x{Ta)/Tcr, is relatively insensitive to ms for fixed 02, 
but increases monotonically with decreasing 64. This situa- 
tion is analogous to the pure SM case where baryon number 
preservation becomes more effective for smaller Higgs quar- 
tic self-coupling. Unlike the SM, however, the YP mass in the 
H phase does not depend on its quartic self-coupling, so one 
may vary at will without encountering phenomenological 
constraints associated with ms. 

The value of my, does, however, affect the strength of the 
second step. Importantly, we find that wherever the two-step 
phase transition occurs in the parameter space, the second step 
of the phase transition is always strongly first order and com- 
fortably satisfies the required BNPC in the Higgs direction. 
As one moves towards lower values of holding 02 and 64 
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FIG. 4: Extremum points of the potential as a function of temperature 
for two choices of model parameters. Curves are labeled according 
to phases defined in Fig. [T] Upper panel: the there is only one crit- 
ical point, corresponding to a SM-like phase transition (O — > H); 
lower panel: the system exhibits two critical temperatures, favorable 
for baiyogenesis, corresponding to transitions at critical temperatures 
T,: O^EandT^: E ^ ff. 



fixed (approaching the left edge of the blue region in Fig. |5]l, 
the order parameter v{Th)/Th of the second step of the phase 
transition dramatically rises. However, this raises the issue of 
whether substantial supercooling of the S-phase would occur 
leading to significant entropy injection, effectively diluting the 
baryon asymmetry through reheating; or if the second phase 
transition ever completes at all. 

To estimate the extent to which supercooling occurs, we 
carry out a numerical calculation of the nucleation rate in the 
semi-classical approximation. We use a modified version of 
CosmoTransitions software ll22ll which rapidly solves for the 
'critical bubble' in theories with multiple field directions'. In 
spirit of maintaining gauge-independence, the code was mod- 
ified to analyze just the 0{T'^) part of the high-T expansion 
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FIG. 5: Phase transition order parameters: Red contours indicate 
constant x{Ta)/Tcr for first step. For the second step, dashed blue 
contours correspond to v{Th)/Th with values {1.2, 1.6, 2.0, 2.4, 
2.8, 3.2, 3.6, 4.0}, read right to left. Outside the contoured region, 
the EWPT proceeds in a single SM-like step (O — >■ H) and is unfa- 
vorable for baryogenesis. 
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FIG. 6: Critical temperature (black), nucleation temperature (red) 
and entropy dilution factor As/sew as a function of triplet scalar 
mass ms for the second step of the phase transition (E H). For 
light triplet masses, the phase transition never occurs (no nucleation 
temperature), as the universe is locked into the meta-stable E-phase 
at zero temperature. Hence, the curves for Tat and entropy dilution 
terminate. 




of the thermal potential (|23] i: 



' A critical bubble is one that continues to expand after formation. 



64(7' 



1 



4 



(25) 
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Here the coefficients 

Da 



mf/ = 125 GeV, ^4 independent 



564 



'iff + ^Vt 
+ ^91) 



202) 



are all gauge-independent. Once the nucleation temperature is 
found, we derived the entropy density injection from the free 
energy: 



As = - ^ 



uv 




dV 






s 


dT 


J 



(26) 



normalized by the entropy density at the time of the elec- 
troweak phase transition 



Sew 



27r^ 
45 



(27) 



where we take = 100. Numerical results are shown in 
Fig. |6l for two choices of 64, and fixed mn and a2. The up- 
per sections of each plot displays the critical temperature 
(black curve) and the nucleation temperature T/v (red curve) 
as a function of the triplet mass ms. As the triplet mass riis 
is parametrically lowered from high values, the S phase be- 
comes more supercooled and the onset of the second phase 
transition is delayed. Below a critical value for m^, the sec- 
ond step of the phase transition never occurs, marked by the 
end of the red curve; the system remains in the metastable S 
phase at zero temperature^. 

In the lower section of each plot, the entropy dilution factor 
As/sew is shown. Entropy dilution typically does not exceed 
3%, which is lies within the typical limits of numerical uncer- 
tainty for baryogenesis computations. 



VI. CONNECTION TO COLLIDER PHENOMENOLOGY 

Because the dynamics of the two-step transition depend on 
the interaction between the isospin triplet and Higgs doublet, 
it is interesting to ask how the viability of this scenario is af- 
fected by the recent observation of a bosonic resonance near 
125 GeV lEl m that may be identified by the SM Higgs. 
Of particular relevance is the report of a marginally signif- 
icant excess izsll of if — 77 decay channel over the other 
expected channels predicted by the Standard Model. As dis- 
cussed in iU, the coupling of the Higgs doublet to the real 
triplet through 02 modifies the branching fraction to two pho- 
tons since the charged components contribute to the rate at 



' Note that a critical temperature T/j may still exist, indicating degeneracy of 
the two minima, even if a transition from one to the other does not actually 
occur. 
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FIG. 7: Predictions for the — 77 branching fraction shift S as 
defined in l |30l ), as a function of triplet mass ms for select values of 
H-E coupling 02. There is enhancement fornegative 02 and suppres- 
sion for positive a2 . 



the one-loop level. In this section we discuss the connection 
of the 02 parameter to the electroweak phase transition. 

The H 'yy branching fraction in the Standard Model is 
dominated by contributions from the W gauge boson and top 
quark loops 12^: 



1 1 

47r niH 



gwFi 



9tFi/2 



(28) 



The W boson gives a positive real contribution, whereas the 
top quark gives negative real contribution. Near the observed 
Higgs mass of ttih ~ 125 GeV, the two processes interfere 
destructively, with the W boson contribution dominating that 
of the top quark. 

The presence of the charged triplet in this model 
modifies iZTll the decay formula by the addition of the scalar 
contribution to the amplitude in i 



pS SM 



1 1 

An rtiH 



where gs 



+ 9tFi/2 ( 
s2woa2/(47r2) 



.9E^o(- 



(29) 



is the effective coupling, inde- 
pendent of 64. The loop function Fq{x) is real and negative 
for ms below threshold x < 4, leading to an additional de- 
structive interference against the W contribution for positive 
a2, and constructive interference for negative 02. We define 
the if — >^ 77 branching fraction shift relative to the Standard 
Model prediction 



5 = 



pS SM _ pSM 
pSM 



(30) 
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FIG. 8: Contours of constant 5 defined in iSOl . superimposed on tree- 
level vacuum stability plot. The yellow and blue regions are color 
coded in the same way as in Fig. |2l and the hashed red area is the 
region where a two-step phase transition favorable for baryogenesis 
is expected. Note the strong correlation between the sign of 02, 5. 

and illustrate the effect of the presence of the charged scalar 
for specific values of 02 in Fig. [T] Reports from ATLAS and 
CMS both mildly favor S > 0. 

In Fig. [8] for fixed 64 and mn, we show the regions in 
parameter space satisfying tree-level vacuum stability with 
the same blue and yellow color-coding used in Fig. |2] The 
red hashed area corresponds to the region in parameter space 
where the desired two-step phase transition is expected to oc- 
cur. We then superimpose contours of constant 5 of ( l30b . Two 
key points emerge: (1) regions where the two-photon branch- 
ing fraction is enhanced occurs for negative 02, where the EW 
vacuum is unstable already at tree-level. (2) Regions where 
a two-step phase transition is likely to occur coincides with a 
reduction of the two-photon branching fraction relative to the 
Standard Model. However, the size of this reduction can be 
lessened for larger triplet scalar mass ms within the region 
where a two-step phase transition occurs. 

Alleviating the tension with the present r{H — > 77) results 
could be achieved in a number of ways. The presence of ad- 
ditional charged degrees of freedom, such as new fermions, 
could compensate for the loop contributions. Alternately, 
one might imagine a two-step transition arising for 02 < 0. 
In the present minimal extension, doing so is not achievable 
since it would imply a tachyonic mass for S° and [see 
©J. But by suitably coupling the S to a new scalar singlet S 
through the operator ^028'^ ■ E, while allowing 5* to obtain 
a zero-temperature vev xs, the new mass relation for S° 

'^1 = -f^s + 7;o.2Vo + cxx% (31) 



could yield m| > for 02 < and C2 > 0. The viability 
of the two-step transition in this case will be the subject of 
forthcoming work. 

VII. DISCUSSION AND CONCLUSIONS 

The observation of a potentially fundamental scalar at the 
LHC makes the paradigm of scalar field-driven symmetry 
breaking in the early universe more realistic than ever. It 
is then interesting to explore possible patterns of symmetry- 
breaking that could arise in the presence of additional scalar 
fields, particularly if these fields and the associated phase tran- 
sition dynamics can help account for the origin of the visible 
and dark matter of the universe. In this study, we have ana- 
lyzed a simple SM scalar sector extension that gives a proto- 
type for a multi-step EWPT wherein the baryon asymmetry 
may be generated prior to the final transition to the SM EW 
vacuum while yielding a contribution to the dark matter relic 
density. This Z2SSM does not provide a complete solution to 
the baryogenesis and dark matter problems, as it lacks addi- 
tional sources of CP-violation required for baryogenesis and 
since for a sub-TeV mass for the new scalar, only a fraction of 
the dark matter relic density can be achieved. Nonetheless, it 
illustrates some of the generic features that may be present in 
more complete scenarios of this type: 

• A strong first EWPT during an initial step in which bub- 
ble nucleation can occur 

• Quenching of the sphaleron transitions during this step 
that results from the coupling of the new scalars to the 
gauge sector of the SM and that is relatively insensitive 
to the mass of the new scalar 

• Subsequent transition to the EW vacuum that does not 
re-activate the sphalerons or lead to dangerous entropy 
injection and whose character is governed by the "Higgs 
portal" scalar operator 

• The presence of a discrete symmetry that ensures stabil- 
ity of the neutral component of the new scalar multiplet 
in the EW phase, making it a contributor to the dark 
matter relic density. 

Many, if not all, of these features can be realized in SM exten- 
sions with higher-dimensional scalar representations, some of 
which may contain new CP-violating phases in the scalar po- 
tential (or through couplings to fermions) as needed for baryo- 
genesis and/or additional dark matter candidates. 

From a phenomenological perspective, it is also an inter- 
esting time to explore the possibilities for realizing this sce- 
nario. As we found in this work, the indications of an excess 
of Higgs decays to two photon final states compared to SM ex- 
pectations creates some tension with the viability of our pro- 
totype scenario, though this tension may be alleviated through 
introduction of new fields. More generally, the presence of 
additional scalar states would lead to novel collider signatures 
that may be explored with increased luminosity and higher 
energy at the LHC. 
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Appendix A: Explicit form of the effective potential 



The zero temperature (Coleman-Weinberg) and finite temperature parts of the 1-loop effective potentials used in the phase 
transition analysis is given here. The dependence on the gauge parameter ^ is explicitly shown. 
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where "free" represents a free-field subtraction. 

Here, the field dependent scalar masses are given by: 

■m%{h,a) 
mQ{h, a) + £,rn^{h, a) 



-^^2 + A/i2 




2 ^ 



-Ms 



2 7,2 I z 

g h + g <7 



2_2 



and the eigenvalues of 

m\{h,a) + ^m'^{h,a) = 
The field dependent gauge boson masses are given by 

ml{h,a)^^ig'+g'')h^ 
m'^{h, cr) = 

and among the SM fermions, we include only the top-quark, whose field-dependent mass given by 



4a - y^ha 



1 



(A3) 
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